Europhysics Letters 



PREPRINT 



Point force manipulation and activated dynamics 
of polymers adsorbed on structured substrates 

P. KRAIKIVSKI, R. LlPOWSKY, AND J. KlERFELD 

Max- Planck- Institut fur Kolloid- und Grenzfldchenforschung, 144-24 Potsdam, Germany 



PACS. 82.35.Gh - Polymers on surfaces; adhesion. 

PACS. 87.15.-v - Biomolecules: structure and physical properties. 

PACS. 87. 15. He - Dynamics and conformational changes. 



Abstract. - We study the activated motion of adsorbed polymers which are driven over 
a structured substrate by a localized point force. Our theory applies to experiments with 
single polymers using, for example, tips of scanning force microscopes to drag the polymer. We 
consider both flexible and semiflexible polymers, and the lateral surface structure is represented 
by double-well or periodic potentials. The dynamics is governed by kink-like excitations for 
which we calculate shapes, energies, and critical point forces. Thermally activated motion 
proceeds by the nucleation of a kink-antikink pair at the point where the force is applied 
and subsequent diffusive separation of kink and antikink. In the stationary state of the driven 
polymer, the collective kink dynamics can be described by an one- dimensional symmetric simple 
exclusion process. 



Introduction. - The thermally activated escape over potential barriers under the influence 
of an external force has been first solved by Kramers for a point particle [1]. Since then this 
process has been extensively studied not only for point particles [2] but also for extended ob- 
jects such as elastic strings. Examples are provided by condensed matter systems: dislocation 
motion in crystals [3,4], motion of flux lines in type-II superconductors [5], or charge-density 
waves [6]. An analogous problem is the activated motion of polymers over a potential barrier, 
which has been considered both for flexible [7] and semiflexible polymers [8,9]. 

In all of these previous studies, the thermally activated motion is induced by spatially 
uniform forces which are applied to the whole polymer or elastic line. In contrast, in the 
present article, we will address the thermally activated motion of polymers over potential 
barriers in the presence of a point force which acts only locally on the polymer. We will 
consider both flexible and semiflexible polymers. 

Our theoretical study is motivated by experimental advances in the manipulation and 
visualization of single polymers using optical [10] and magnetic [11] tweezers, or scanning 
force microscopy [12]. In ref. [12] it has been demonstrated that these techniques allow to 
experimentally apply localized point forces to a polymer adsorbed on a substrate. Polymers 
that are strongly adsorbed onto crystalline substrates such as graphite or mica experience a 
spatially modulated adsorption potential reflecting the underlying crystal lattice structure and 
giving rise to preferred orientations of the adsorbed polymer. For such systems, the dynamics 
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of the adsorbed polymer is governed by thermal activation over the potential barriers of the 
surface potential. 

One example of polymers adsorbed on a structured surface are self- assembling polymer 
chains consisting of long-chain alkanes and alkylated small molecules on crystalline substrates 
such as the basal plane of graphite [13]. The alkyl chains orient along the substrate axes 
thereby providing an effective periodic adsorption potential. Also biopolymers such as DNA 
or polyelcctrolytcs can be oriented on the basal plane of graphite by using long chain alkanes 
as an oriented template layer [12,14]. It has been demonstrated experimentally that these 
polymers can be manipulated individually on the structured surface by applying point forces 
using the tip of a scanning force microscope [12]. 

Our main results are as follows. At low forces, the dynamics of the polymer is governed 
by thermal activation and nucleation of localized kink-like excitations as shown in fig. ^ We 
calculate the critical point force below which the polymer moves by thermal activation over the 
barriers of the adsorption potential. The steady state of this activated motion determines the 
profile and velocity of the moving polymer and is governed by the (collective) driven motion of 
the kink excitations which can be described as a one-dimensional symmetric simple exclusion 
process of these excitations. Our results for the critical point force, the velocity, and the profile 
of the moving polymer are accessible in manipulation experiments on adsorbed polymers and 
allow to extract material parameters of the polymer and the substrate structure from such 
experiments. We will first present detailed calculations for stiff, semiflexible polymers. Results 
for flexible polymers are discussed in the end. 

Model. - We consider the dynamics of a semiflexible polymer adsorbed to a planar two- 
dimensional structured substrate under the influence of an external point force F p pulling 
the polymer. A generic model of the substrate structure is a double-well potential that is 
translationally invariant in one direction, say the x-axis as in fig-H The semiflexible polymer 
has a bending rigidity k and persistence length L p = 2k /T where T is the temperature in 
energy units. We focus on the regime where the potential wells are sufficiently deep and 
narrow so that the adsorbed polymer is oriented along the ir-axis and can be parameterized 
by displacements z(x) perpendicular to the x-axis with — L/2 < x < L/2, where L is the 
projected length of polymer, see fig. ^ The Hamiltonian of an oriented polymer is given by 



J-L/2 LZ J 

i.e., the sum of bending and potential energy [8,9]. We consider a piecewise harmonic double- 
well potential 



with Vq(z) = t;Vq(\z\ — a) 2 , where Vo is the depth of the potential. The potential J2J contains 
the action of a point force pulling the polymer at the point x — x p with a force F p in the 
z-direction. For zero point force F p = 0, the potential is symmetric, translationally invariant 
in the x-direction, has a barrier height Vb& 2 /2, and the distance between minima is 2a. For 
F p > 0, the point force in J2J breaks the translational invariance of the system. 

Our assumption of an oriented polymer is valid if U-turns of the polymer within a single 
potential well are suppressed by the bending energy. This is the case if the size 2a of each 
potential well in the z-direction is smaller than the persistence length L p . This condition 
is typically fulfilled for adsorbing substrates structured on the nm-scale [13]. Furthermore, 
the polymer should be strongly adsorbed, which corresponds to a small density of thermally 
induced kink excitations, i.e., Ey. ^> T where Ef. is the kink energy, see eq. (QJ below and 




(1) 



V p (x, z) = Vq(z) - F p S(x - x p )z 



(2) 



ref. [8] . 
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Fig. 1 - (a) Kink-antikink configuration of a semiflexible polymer in a double-well potential V 
under the action of a point force F p displacing the midpoint in the ^-direction to a value z m . The 
configuration Zk(x) is calculated for F p /F c = 0.19, z m /a — 0.21 (Li/wk = 1-1, Li/wk = 15) and has 
an energy E/2Ek = 0.72. (b) The midpoint z m (in units of a) as a function of the external force (in 
units of critical force 4Ek/a). (c) Energy E(z m ) (in units of 2Ek) of a kink-antikink pair as a function 
of the midpoint z m (in units of a) for different forces F v /F c — 0, 0.19, 0.5. 



The overdamped motion of the polymer is described by [8, 9] 

yd t z = —— + C(x,t) = -Kdiz-V£(z) + Fp5(x-x p )+((x,t) , (3) 

where 7 is the damping constant and £(x, t) is a Gaussian distributed thermal random force 
with (C) = and correlations (((x, t)((x' , t')) — 2^T5{x — x')8(t — t'). We neglect longitudinal 
motion of polymer segments (sec rcf. [9] for a discussion) and do not study the effects of an 
external tension or compression. For V p — 0, tension and compression have been considered 
in [15] and [16], respectively. 

Static kinks. - First, we calculate the stationary shape of the semiflexible polymer that 
is deformed by a point force acting at its midpoint into a kink-antikink configuration Zk(x) 
as shown in fig.^Ji. This configuration is obtained by displacing the polymer at the midpoint 
where the point force acts to a prescribed position z m and letting the rest of the polymer 
equilibrate. Therefore, we have to solve the saddle-point equation 8H/Sz — for the energy 
i.e., eq. for the time-independent case and in the absence of noise (C = 0), with 
appropriate boundary conditions and a prescribed position Zk{x p ) — z m . For z m > the 
kink configuration crosses the barrier at two points, see fig. we choose the origin x = 
and the length L2 such that these points are Zfc(0) = and z^iE^) — 0. The polymer has 
a total length L = L\ + L2 and extends from x — —Li/2 to x = Li/2 + L2, and the force 
acts at the midpoint x p — L2/2. The kink-like configuration has to fulfill four boundary 
conditions, z k (~L 1 /2) = z k (-Li/2 + L 2 ) = -a and z' k \- Ll /2 = z k \^ Ll / 2 +L 2 = 0- At the 
midpoint x p — L2/2, we fix the displacement z m of the polymer Zk(x p ) = z m , and the point 
force causes a discontinuity in the third derivative, z' k "(x p +) — z' k "(x p —) — F p /k. In addition, 
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Zk{x) and its first two derivatives have to be continuous at the midpoint, and Zk(x) and its 
first three derivatives have to be continuous at each crossing point xq = 0, £2- 

Away from the point force, i.e., for x 7^ x p the saddle point solutions are linear combi- 
nations of the four functions exp(±x/wk) exp(±ix/wk) where Wk = V^in/Vo) 1 ^ 4 is the kink 
width. Construction of the solution through the four regions separated by the crossing points 
and the midpoint then requires to determine 16 linear expansion coefficients and the two pa- 
rameters Li and z m as a function of the system size L and the remaining model parameters 
including the point force from the boundary and matching conditions. The resulting shapes of 
the kink-like polymer configurations are shown in fig.^i. fig. ^ shows the energies E(z m ) of 
the kink- like configuration as a function of z m for different point forces F p . For low forces the 
energies E(z m ) in fig.^ have two stationary points, a stable minimum at z m = z mmin < 
(the midpoint does not cross the barrier) and an unstable maximum at z m — z rn , nuc > 0. This 
maximum is unstable with respect to further displacement of the midpoint and represents the 
critical nucleus configuration. For F p = 0, we obtain another stable minimum at z m = a (the 
midpoint reaches the next potential well) which is the static kink-antikink solution [8]. The 
width Wk of a static kink and its characteristic energy E)~ are given by 

™ fc = V2( K /F ) 1/4 , E k = a 2 K^V 3/4 /V2 . (4) 

In the limit of large L, we can find analytic expressions for the resulting stationary positions 
z m ,min and z m<n ,uc as a function of the applied force F p , see fig.^D. We find that there are no 
stationary positions if the point force F p exceeds a critical value F c given by 

F c = AE k /a = 2V2clk 1/4 Vo /4 . (5) 

The midpoint displacement z m ^ m i n < in the stationary minimum is a linear function of the 
external force, z m ,mm = —a (1 — F p /F c ) and reaches the barrier at z m , m in — for F p — F c , see 
fig. \\]p. This force-displacement relation describes the linear response of the polymer before 
crossing the barrier. For the midpoint displacement in the unstable nucleus configuration 
z m nuc > 0, on the other hand, we obtain the following set of two equations for z m nuc and 
L2, 

F p /F c = (cos x - sin %)e~ x \ x=L2 / 2Wh , z m ,rmc/u = 1 - (sina; + co$x)e~ x \ x=L , 2w ^ . (6) 

As shown in fig.^D, F p is decreasing for increasing z m ,nuc as the critical nucleus configuration 
widens for small point forces. The negative values of F p for large z m%nuc indicate that for a 
semiflexible polymer the kink-antikink configuration reached for z rn = a is stabilized by an 
energy barrier. Only below a negative threshold force F~ = —F c e~^l 2 < the kink-antikink 
configuration becomes unstable. 

Kink nucleation. - Now we turn to the activated kink nucleation in the presence of a point 
force pushing the polymer over the potential barrier. The point force breaks the translational 
invariance in cc-direction and kink-antikink pairs are only nucleated at x = x p with a rate J per 
unit time. This thermally activated process is governed by an energy barrier which is given 
by the excess energy AE n of the critical nucleus configuration. The energy of the critical 
nucleus can be obtained from the energy profiles E(z m ) shown in fig. ^ as the difference 
AE n = E(z m ,nuc) — E{z m ^nin) between minimum and maximum values of the energy E(z m ) 
of the kink-like configuration. We find AE n ~ 2E k (1 — F p /F c ) 2 , which vanishes as the force 
approaches the critical value F c . The activation energy enters the nucleation current 

J = {Qn/2ir) exp (-AE n /T) with Q\ = |u; n ,o|w s ,o JJ(w s ,p/w„,p) , (7) 

p>0 
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Fig. 2 - (Left) The shape of a semiflexible polymer pulled over a periodically structured surface by 
a point force acting at the midpoint. The horizontal lines indicate the position of potential barriers. 
The thick solid line shows a typical polymer configuration z(x), the thin line the average shape (z(x)). 
(Right) The stationary kink density pt (x) as a function of the distance from the point x v where the 
force is acting on the polymer. 



which shows Arrhenius-type behaviour. The prefactor Q n includes the spectrum of attempt 
frequencies L/J n ,p and uj StP (p — 0, 1, ...) for phononic fluctuations around the critical nucleus 
configuration and the straight configuration z m — —a, respectively. We find one unstable 
negative mode w n . < 0, which diverges as w„ i0 = (Vq/7)[1 — 2 4 / 3 (l — F p /f CiK ) -4 ' 3 ] upon 
approaching the critical force F p ss F c , a bound state with < u> n i < V0/7, and a set 
of positive modes u> n>p > Vq/j with the same level spacing as the modes of the straight 
configuration. It is important to note that two translational modes (for kink and antikink) 
only exist if the point force is zero because the point force breaks the translation invariance. 
Close to the critical force F p ~ F c , we obtain Q 2 n « (Vb/7) 2 [1 - 2 4 / 3 (l - Fp/F c )" 4 / 3 ]. 

Collective kink dynamics. - After nucleation of a kink-antikink pair at x = x p by thermal 
activation, kink and antikink are driven apart by a small force ~ F\ z /wke~ L ' 2 l Wk , which decays 
exponentially with the distance Li > Wk between kink and antikink. This exponential decay 
is characteristic for a point driving force which interacts only over a distance ~ Wk with the 
kink and very different from the case of a spatially uniform force, where kinks experience a 
spatially uniform driving force [8,17]. For separations Li > Wk the kink diffuses essentially 
freely with a diffusion constant Dk — ITwkjZ^a? [8]. 

A spatially localized driving force also leads to a distinct steady state motion of the polymer 
in a periodically continued potential, see fig. |3 This motion can be described in terms of the 
collective dynamics of an ensemble of kinks and antikinks which are generated at the single 
point x — x p by the point force and subsequently separated by the exponentially decaying 
force. For the following discussion we choose coordinates such that x p = 0, and the polymer 
extends from —L/2 < x < L/2. Because a point force creates kink-antikink pairs only at 
x = 0, we find an ensemble consisting only of kinks in the region x > and an ensemble 
consisting only of antikinks in x < 0. As two (anti-)kinks have a mutual short-range repulsion 
of range Wk, we have an ensemble of diffusing kinks (antikinks) with a hard-core repulsion on 
the interval L/2 > x > (—L/2 < x < 0). In order to treat the non-equilibrium dynamics 
of these ensembles, we introduce a discrete one-dimensional lattice of possible kink positions 
with spacing Ax = Wk which allows to map the dynamics of each ensemble onto the symmetric 
simple exclusion process (SSEP) with open boundaries [18,19]. In the following we consider the 
kink ensemble (x > 0); the antikink ensemble (x < 0) can be treated analogously. In the kink 
ensemble, the kink particles are freely diffusing, i.e., they have symmetric rates D = D^/w\ 
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for hopping to the right and left on the lattice X{ = iw k [i = l,...,iV with N — L/2w k ); they 
interact through their hard-core repulsion. In the SSEP, boundary conditions are specified by 
rates a and 6 for particles to enter the system at the left (i = 1) and right (i = N), respectively, 
if that site is empty. For the kink ensemble we have a — J, as kinks are nucleated at i = 1 
with the Kramers rate JJJ, and S — as no kinks enter the system at i = N. Furthermore, 
kinks leave the system diffusively, at i — 1 by annihilation with an antikink and at i = N by 
relaxation of the free polymer end. 

Despite the hard-core interaction the stationary density profile p k {x) of kinks in the SSEP 
fulfills the stationary diffusion equation, d 2 p k — [18,19]. Furthermore, our boundary condi- 
tions are equivalent to boundary conditions p k (0) = w^ 1 min(a/-D, 1) and p k (L/2) = for the 
stationary kink density at the ends of the system. For a > D the system reaches its maximal 
kink density w^ 1 at x = 0. The resulting linear density profile p k (x) is 

Pk(x) = p fe (0)(l - 2\x\/L) with p fc (0) = w- 1 wm(a/D, 1) = min( Jw k /D k , l/w k ) (8) 

as shown in fig. [21 (right). The average distance between kinks is l/p k (x) and at each kink the 
polymer position changes by Az = —2a leading to a characteristic parabolic polymer shape 
(z(x)) — z m = —2a dxpk(x) = — 2(a/wk) mm^Jw 2 ,/ Df., — \x\/L) in the stationary 

state as shown in fig. [21 (left). The average velocity v z = (dtz) of the polymer in the z- 
direction is determined by the stationary current Jssep = —D^dxpk — min(J, Dk/wfywk/ L 
of the SSEP. Only for small nucleation rates J <C Dk/w\ the kink interaction can be ne- 
glected and the current is directly given by the Kramers rate {JJ, Jssep ~ Jwu/L. Dur- 
ing the time 1/ Jssep the polymer advances by a distance 2a leading to v z = 2a Jssep ~ 
2amin(J, D k /w\)wk/ L. 

Flexible polymers. So far we considered semiflexible polymers dominated by their 
bending energy. In this section we want to outline the main results for flexible Gaussian 
polymers governed by entropic elasticity with a tension a — 2T/b where b is the Kuhn length. 
The Hamiltonian of a flexible Gaussian polymer on a planar two-dimensional substrate is 
given by 

f L c /2 

H= ds 

J-L a /2 

where we integrate over the arc length s with —L c /2 < s < L c , and L c is the contour length 
of the polymer. The translationally invariant potential V{z) is a function of z only. Therefore 
fluctuations in the x-coordinate decouple and are Gaussian with moments ((x(L c ) — x(0))) = 
and ((x(L c ) — a;(0)) 2 ) w L c b/2. The Rouse dynamics of the z-coordinate of the polymer is 
given by 

7 9 t z = ad 2 s z - V£(z) + F p S(s - s p ) + C(s, t) , (10) 

where 7 is the damping constant and £(x, t) is a Gaussian distributed thermal random 
force. The point force is on the monomer s = s p . For a flexible polymer the kink width 
is w>fc,er = (c/Vo) 1 / 2 and the kink energy E k ^ = a 2 (aVo) 1 ^ 2 [7,17]. As for the semiflexible 
polymer we can calculate the energy E(z m ) of a kink-antikink configuration with prescribed 
midpoint z k {s p ) — z rn . For a flexible polymer the displacements in the stationary minimum 
at z m = z m _ m i n and the maximum representing the nucleus with z m — z m ^ nuc are both lin- 
ear functions of the external force, z m ,min = —a (1 — F p /F c ^) and z m ,nuc = a(l — F p /F c ^), 
where the critical force for the flexible polymer is given by F Ci<T = 2E kyCr /a — 2a(aVo) 1 ^ 2 ■ 
The nucleation current J CT for the flexible polymer is given by the same expression as for a 
semiflexible polymer with the excess energy AE n ~ 2E ka (1 — F p jF ca ) . The spectrum of at- 
tempt frequencies uj n p for the critical nucleus shows slightly different behaviour for the flexible 



^[{d s x) 2 + (d s z) 2 ]+V{z)\ , (9) 
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polymer as the unstable negative mode Lo n ,o ~ — 3Vb/7 does not diverge for F p « F Cj(J , and we 
finally obtain Q n w v / 3Vb/7. The collective kink dynamics for a flexible polymer can also be 
mapped onto a one-dimensional SSEP. As a function of the arc length s, we find a linear sta- 
tionary kink density profile pk(s) = /9fc(0)(l — 2\s\/L c ) with pk(0) = min(Jo-u;fc,o-/-Dfc,o-, l/wk, a ) 
and a parabolic shape (z(s)) = —2(a/wk,o-) rnin( J a w%. a /Dk,o-, — |s|/L c ) analogously 

to the semiflexible polymer, cf. eq. (JSJ, where -Dfe.o- = Fwk,a/lo? [17] is the kink diffusion 
constant of the flexible polymer. In the real space coordinates of the substrate, however, 
the resulting shape is ((x(s)} , (z(s))) = (0, (z(s))) and thus, the parabolic shape is lost 
due to the decoupled Gaussian fluctuations in the x-direction. The result for the velocity 
v z = 2aJsSEP ~ 2amin(J CT , D). >a /vyt a )wk.a / L c is analogous to the semiflexible polymer. 

Conclusion. - In summary, we described the activated motion of single adsorbed poly- 
mers on a structured substrate displaced by localized point forces, which can be realized 
experimentally using, e.g., scanning force microscopy tips. The dynamics is governed by kink- 
like excitations for which we have calculated shapes, energies, and critical point forces. Kink 
and antikink pairs are locally nucleated by the point force and then undergo a separation 
which is diffusive on separations larger than the kink width Wk- We have calculated the nu- 
cleation rate (JJJ using Kramers theory. The collective kink dynamics can be mapped onto a 
one-dimensional symmetric simple exclusion process (SSEP). Using this mapping we find the 
average polymer velocity and a characteristic average parabolic shape for a driven semiflexible 
polymer. 

* # * 

We thank S. Klumpp for discussions on the simple symmetric exclusion process. 
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